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Abstract 

We demonstrate numerically that proton-proton (pp) scattering observables can be determined directly 
by standard short range methods using a screened pp Coulomb force without renormalization. In examples 
the appropriate screening radii are given. We also numerically investigate solutions of the 3-dimensional 
Lippmann-Schwinger (LS) equation for a screened Coulomb potential alone in the limit of large screening 
radii and confirm analytically predicted properties for off-shell, half-shell and on-shcll Coulomb t-matrices. 
These 3-dimensional solutions will form a basis for a novel approach to include the pp Coulomb interaction 
into the 3N Faddeev framework. 

PACS numbers: 21.45.-v, 21.45.Bc, 25.10.+S, 25.40. Cm 
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I. INTRODUCTION 



The action of the Coulomb force in pp scattering can be rigorously treated using the Vincent- 
Phatak method l|. We propose an alternative manner using a screened Coulomb force, despite the 
well know fact that the screening limit does not exist. Namely the pp on-shell scattering amplitude 
acquires an oscillating phase factor if the screening radius goes to infinity 2], 3, 4, 5|. This phase 
factor is known and can be removed, a step known in that context under the name renormalization. 
However, as we shall show, if one is interested in the pp observables (not in the phase shifts) where 
that phase factor drops out, all scattering observables can be obtained in the standard framework 
of short range interactions. This will be demonstrated in section II for suitably chosen screening 
radii. 

In view of a forthcoming paper [(J related to the pd system, we further investigate in section 
III properties of the screened 3-dimensional Coulomb t-matrix < p '\t^(E)\p >. This t-matrix 
is a solution of the 3-dimensional 2-body LS equation driven by the screened Coulomb potential. 
Namely to catch the full action of the Coulomb force in the pd system a partial wave truncated pp t- 
matrix is insufficient and the complete 3-dimensional Coulomb t-matrix has to be used. Analytical 
properties of that screened Coulomb t-matrix, off-the-energy-shell, half-shell and on-shell have been 
studied in the past . These investigations, however, mostly rely on insights gained 

for fixed partial wave states. The mathematical rigor in the summation of the partial wave sum 
to infinity leaves room for improvement. Therefore we felt that a numerical study is justified to 
verify statements given there: the screening limit of < p '\t^(E)\p > exists for ^— ^ E ^ and 
coincides with the unscreened pure Coulomb force expression, which is known analytically 



and references therein; that screening limit exhibits a discontinuity if p approaches y / rripE, E > 
from above or below; the screening limit of the on-shell t-matrix < pp'\t^(E = ^r)\p > approaches 
the analytically known unscreened Coulomb on-shell t-matrix up to a given infinitely oscillating 
phase factor. Here we want to numerically investigate at which R-values these limits are reached 
with adequate accuracy. We conclude in section IV. 



II. THE ON-SHELL PP T-MATRIX WITH SCREENED COULOMB POTENTIAL AND 
THE PP OBSERVABLES 

Let be the screened Coulomb potential between 2 protons normalised such that turns 
into the pure pp Coulomb potential for R, the screening radius, going to infinity. Together with 
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the strong interaction V this determines the 2-body pp t-matrix via the LS equation 

t = V + V R + (V + V R )G t , (1) 

2 

where Go is the free propagator. That equation is solved at the pp cm. energy E = ^- projected 
on a set of partial wave basis states \p(ls)jm;tmt >, with p,l,s,j and m the relative momentum, 
orbital angular momentum, total spin, total angular momentum and its magnetic quantum number. 

The total isospin quantum numbers for two protons are t = 1 and mj = —1. This leads to the 
on-the-energy-shell t-matrix element 

< p(l's')j'm'\t\p(ls)jm >= 5 s > s 5 j > j 5 m i m t s l ? l (p,p) , (2) 

where the Pauli principle dictates (— ) l+s = 1 and we took s to be conserved. 
The full 3-dimensional antisymmetrized on-shell t-matrix is given as 

< p 'm^m^l^l — P\2)\pm\'m2 > , (3) 

where mi (m! i ) are the individual spin magnetic quantum numbers and p = pp, p' = pp ' the initial 
and final relative momenta. 

The standard partial wave decomposition leads to 

< p 'm'im'2\t(l — Pyi)\pm\m<i >= ^^(^ 2 S ' m 'i m 2 m s)(2 2 s ' m i m 2 m s) 

s 

00 j j+s j+s 

EE E E E( //s ^' m - m ) y ^ ')$(M(i + (-) l+s ) 

j=0 m=—j V=\j-s\ l=\j—s\ m[ 

^2(lsj,mim s ,m)Y l l ll (p) . (4) 

The strong force can be neglected beyond a certain j max and there only the screened Coulomb 
t-matrix t R is present, which is diagonal in I and independent of s and j. In a well known manner 
one adds and subtracts a finite sum up to j ma x with t R only and this completes the infinite sum over 
j containing only t R . That infinite sum is identical to the 3-dimensional antisymmetric screened 
Coulomb t-matrix. Thus turns into 

<p'm' l m' 2 \t{l - P 12 )\pm 1 m 2 >= S^J^^ <p'\t R \p > -5 m / m2 <5 m > 2mi < p '\t R \ -p> 

^-^ 00 j j+s j+s 

+ ^(2 2 s,m ' im ' 2m ' s ^2 2 s,mim2ms ^^2 E E E E( //sj ' m '' m 's' m ) 

s ' j=0 m=-j l'=\j-s\ l=\j-s\ mj 

Y Vm [{P ')(4(P,P) - <Wd)(l + (-) /+s )^(/5J,m / m s ,m)y4 ; (p) . (5) 

mi 
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Now as is well known 



BE 



4J the limit of that expression does not exist for R — > oo. In that limit 
each term in ([5]) acquires the same infinitely oscillating factor e 2 **^), where &r(p) is given below. 
If one is interested in scattering phase shifts it is unavoidable to keep track of this oscillating factor 
which in that context runs under the name renormalization j|. However, if one is interested in the 
pp observables, the cross section and all sorts of spin observables (note &r(p) is independent of 
spin magnetic quantum numbers), where the on-shell t-matrix appears together with its complex 
conjugate, the oscillating factor drops out. In that case one does not even has to know the analytical 
form of <&r(p). It is sufficient to know that the limit of large screening radius generates just a phase 
factor. 

This is the main message of this section: the pp observables based on the strong and the 
screened Coulomb force can be calculated without renormalization using standard short range 
methods. Though not explicitly stated in [j], J] this insight is in the spirit of these authors. It 
remains to establish the values of the parameter R at which the observables get independent of R. 

Based on ([5]) all pp-scattering observables are given by well known analytic expressions [lo| . 
here, however, we use a more modern nomenclature for the various spin observables [ill ]. We use 
the following screening form which depends on two parameters, the screening radius R and the 
power n: 



V c \r) = e -e-^ . (6) 



At a given value n the pure Coulomb potential results for R — > oo. We use n = 1, 2, 3, and 4. 
As has been shown in [2] based on 3, 4], the related phase <&r(p) is given as 

$r(p) = -r}[ln{2pR) - 7 /n] (7) 

2 

where 7 = 0.5772 ... is the Euler number and r\ = the Sommerfeld parameter. 
Considering only the screened Coulomb force ([6]) the leading term in ([TJ is given by 

e 2 1 

<p'\V^\p>= — ^7^ — ^77 / drsin(\p- p'\r)e~^ n . (8) 
2vH \p-p'\ J 

By a simple partial integration it can be shown that this matrix element approaches the n- 
independent limit 

e 2 

lim <p'\V c R \p>= . (9) 

R^oo 2n 2 \p -p'\ 2 

Therefore for suitably large R- values the leading term and thus the solution of the LS equation (JTJ) 
will approach a n-independent limit. 



In Figs. Q] and [2] we demonstrate at E l ^ b = 13 MeV for several pp observables independence on 
n for a suitably large R-value and the perfect agreement to the exact Vincent-Phatak results at 
E l p b = 13 MeV. The deviation with respect to the standard Vincent-Phatak approach [l|, which 
treats the pp Coulomb force rigorously, for different pp observables and values of n = 1, 2, 3 and 4 
and screening radius R = 120 fm is under ~ 1%. 

In Figs. [3] and 2] we show the convergence with respect to R for n = 4 for a number of pp 
observables at E l £ b = 13 MeV. The corresponding results for E l ^ b = 50 MeV are shown in Figs. [5] 
and [6J The resulting limiting values agree very well with the Vincent-Phatak results. At E l ^ b = 
13 MeV the limiting .R-value is R = 120 fm and at E l p ab = 50 MeV R = 60 fm. With decreasing 
energy the limiting i?-value increases. 



III. PROPERTIES OF THE 3-DIMENSIONAL SCREENED COULOMB T-MATRIX 



As will be shown in a forthcoming article [6J the 3-dimensional screened pp Coulomb t-matrix 
< p '\tf(E)\p > occurs naturally in a certain type of Faddeev equation of the pd scattering problem. 
There it appears off-the-energy-shell (with the exception of isolated points). For the unscreened 
pure Coulomb force the off-the-energy-shell expression is analytically known 

ELQ. 

We numerically investigate the screening limits of < p'\tf(E)\p > for 2- ^ E ^ 2- (off-shell), 

2 2 

for p' p and E = ^- (half-shell) and for p' = p and E = ^- (on-shell). Here we want to 
numerically investigate at which i?-values the limits are reached with adequate accuracy and how 
they are related to the corresponding unscreened pure Coulomb force expressions. 

To that aim we regard the LS equation for two protons interacting only with the screened 
Coulomb potential V^. The off-shell t-matrix element < p'\t^(E = ^-)\p >= t^(p',P, x = p-p E) 
fulfills for given energy E the equation 

1 f°° I' 1 1 

t?(p',p,x) = —v?(p',p,x,l)+ dfp" 2 \ dx"v?(p',p»,x,x») t?(p»,p,x") (10) 

27r Jo J -i E + ie-2- 

flip 

with 

r-2n 

v^(p',p,x r ,x)= d<pV c R (p',p,x'x + \/l - z'Vl - x 2 cos(t)) . (11) 

J 

Eq. ()10p can be solved after discretizing x and the continuous momentum variables using direct 
matrix inversion or generating the Neumann series and applying Pade summation. 

We solved ([10]) at three energies: E = 3 MeV, 13 MeV, and 50 MeV. For n = 1 the leading 
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term in (jlOp can be calculated analytically 

v?(p',p,x',x) = e " _ . (12) 

7r-i / (p' 2 + p 2 — 2p'px'x + jp) 2 — *ip' 2 p 2 {l — x /2 )(l — X 2 ) 

For n > 1 this is no more possible and a two-dimensional numerical integration is required to 
get the leading term 

e 2 f 2lT f°° 1 _ (11 „ 

v(p ,p, x , x) = — 2 / tf</> / (ir — - e ^ . (13) 



Jo yp 2 + p' 2 — 2pp'(x'xy/l — x l2 \J\ — x 2 cos(f)) 

On the other hand the pure off-shell Coulomb t-matrix is known analytically (0,111 and refer- 
ences therein) 

<^€^>-&w^w (14) 

with 

I{x) = - [ 2 F X (1, mj; 1 + ir,; - 2 F 1 (l, i m 1 + wj; ^j)] (15) 

a; x — 1 x + 1 

a 2 2 2 I"! 

and x 2 = 1 + ^ p ki^2^jpp k ' • 2-^1 is the hypergeometric function [14], which we determined using 
subroutines from [lEf] and paid attention to the vanishing small positive imaginary part of k 2 . 

As an illustration of our general results we show in Fig. [7] the limiting behavior for the real and 
imaginary parts of the off-shell screened Coulomb t-matrix t^(p,p f , x) at E p ab = 13 MeV and fixed p 
and x values as a function of p'. This energy corresponds to the on-shell momentum k = 0.396 fm . 
The small R- values R = 20 fm and R = 60 fm are quite insufficient, especially for the imaginary 
part, to reach the pure Coulomb off-shell values. For the higher i?-values, R = 120 fm and beyond, 
t^(p,p',x) converges very well and the limit coincides, as expected, with the pure Coulomb off- 
shell t-matrix. For p' = k = 0.396 fm -1 one reaches the half-shell point and in its neighborhood 
a discontinuity develops with increasing i?-values, if one approaches k from below or above. That 
discontinuity is very well known to exist for the pure half-shell Coulomb t-matrix and reproduced 
for the convenience of the reader: 



m p 27T Z \p — p | z 1 — e Z7rr ' 



1 k 2 > p' 2 

(16) 



That discontinuity is separately shown in Fig. [8] for both parts of the pure Coulomb t-matrix, 
for predictions based on Eqs [T4TI15I The screened Coulomb t-matrix is, of course, continuous at 
p' = k for each fixed i?-value, but nevertheless the tendency to develop that discontinuity can be 
seen with increasing R. 



Now as has been emphasized in [?], [s| that discontinuity would be absent if the limits p' — > k 
and R — > oo are performed such that |p' — — > 0. To show that numerically a more subtle 
investigation is required, which we did not undertake. 

Fig. [9] exemplifies the situation for the off-shell t-matrix at negative energy E^ b = — 13 MeV, 
where the t-matrix is real. There the limit is reached already around R = 20 fm. 

Let us turn now to the half-shell pure Coulomb t-matrix, which is analytically given by \}\ - 



k2 leu ti'^ Jc^ 1 

< p'\^(ir)\k >- c^^—^-r , (17) 



where q = p ' — k is the momentum transfer, o"o = T(l + irf) is the pure Coulomb phase shift and 
Cq = cxp ^_ 1 is the Coulomb penetrability. The direct comparison of this limit and the screened 
Coulomb half-shell t-matrix is shown in Figs. [10] and [TT] for the real and the imaginary part of t 
at E l p ab = 13 MeV, respectively. On both fi gures, in upper row discrepancy due to the oscillating 
factor e Ri - > is seen. After removing that factor (by procedure called renormalization) the 

screened half-shell t-matrix approaches in the limit R — > oo the half-shell Coulomb t-matrix, what 
is shown in lower row of Figs. [10] and [TT1 The results of renormalization are striking, especially 
for the imaginary part of the half-shell screened t-matrix. The screening radii about R = 60 fm is 
sufficient to describe the pure Coulomb t-matrix, however for smaller angles one has to go even to 
higher i?'s. 

Finally let us regard the screened on-the-energy-shell t-matrix element, which acquires for large 
R- values the oscillating factor e r( - ' . After renormalization with this factor the screened 

on-the-energy-shell t-matrix approaches in the limit R — > oo the Coulomb scattering amplitude 



-energy 



A c (e) 

t c (k,k,x)\ renor malized = e R<K (k,k,x) — > ~— (2^p = m p (2ir) 2 'W 2 sin 2 1 ■ 

This is demonstrated in Fig. [12] for the real part and in Fig. [13] for the imaginary part of the 
on-shell t-matrix at E l p b = 13 MeV. The upper rows show unrenormalized screened t-matrices 
while the lower rows show them after renormalization. It is clearly seen that the renormalization 
is required to get the Coulomb on-shell amplitude (shown by the thick solid line). The effect of 
renormalization is very pronounced for the imaginary part. It is negative without renormalization 
and shows a strong dependence on the screening radii R. After renormalization and for sufficiently 
large R- values it perfectly overlaps with the imaginary part of the pure Coulomb on-shell amplitude. 
At forward angles, which corresponds to x close to 1, R = 180 fm is required to reach good 
agreement with the pure on-shell Coulomb t-matrix. 



IV. SUMMARY AND OUTLOOK 



We numerically solved the 3-dimensional LS equation for a general off-shell screened Coulomb t- 
matrix with different types of screening. That t-matrix taken on-shell together with a finite number 
of partial wave projected t-matrices generated by the sum of the screened Coulomb force and the 
nuclear force and corrected for the partial wave projected screened pure Coulomb t-matrix leads to 
correct pp observables at suitably chosen finite screening radii. The renormalization phases drop 
out automatically in the pp observables since they are products of the full t-matrix and its complex 
conjugate. Thus renormalization is not required and even the knowledge of the renormalization 
phase is not required when calculating observables. 

Finally we numerically checked that the screened 3-dimensional Coulomb t-matrix has the off- 
shell, half-shell and on-shell limits in relation to the general unscreened pure Coulomb t-matrix, 
which have been analytically predicted in the literature. We felt that this numerical 3-dimensional 
investigation supplements well those previous analytical studies, which lack mathematical rigor 
when summing partial wave results to infinite order. 

The resulting 3-dimensional screened general Coulomb t-matrix will be used in a forthcoming 
Faddeev treatment of pd reactions. 
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FIG. 1: (color online) The convergence of the pp scattering cross section (J^), analyzing power (A/), spin 
correlation coefficient (C yy ), and spin transfer coefficient (K% ) at E l p ab = 13 MeV as a function of the cm. 
scattering angle calculated with the screened Coulomb force and the CD Bonn nucleon-nucleon potential 18| , 
which is kept for the partial waves up to j < 3. The screening radius is R — 120 fm and n ~ 1 (dotted line), 
n = 2 (dashed-dotted line), n — 3 (dashed line), and n = 4 (solid line). The curves for n = 1 to n = 4 all 
overlap on the scale of the figure. The exact Vincent-Phatak result is given by thick dots. 
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FIG. 3: (color online) The pp scattering cross section (^), analyzing power (A y ), spin correlation coefficient 
(C yy ), and spin transfer coefficient (Kjj ) at E l ^ lh = 13 MeV as a function of the cm. scattering angle 
calculated with the screened Coulomb force and the CD Bonn nucleon-nucleon potential [18J , which is kept 
for the partial waves up to j < 3. The screened results are for n = 4 and different values of the screening 
radius R: R = 20 fm (dotted line), R — 60 fm (dashed-dotted line), R = 120 fm (dashed line), and 
R = 180 fm (solid line). The thick dots are the Vincent-Phatak's exact results. 
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FIG. 4: (color online) The same as in Fig. [3] for other spin transfer coefficients. 
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FIG. 7: (color online) Limiting behaviour of the real (upper) and the imaginary (lower) parts of the off- 
the-energy-shell screened t-matrix t^(p,p',x) at E l £ h = 13 MeV, p — 0.36 fm _1 , and x — 0.71 as a function 
of the p' momentum for n = 4 and different values of the screening radius R: R — 20 fm (dotted line), 
R = 60 fm (dashed-dotted line), R = 120 fm (dashed line), R = 180 fm (thin solid line), R = 500 fm (thick 
solid line). The pure Coulomb off-shell result of Eq. p^j) is given by thick dots. The half-shell situation is 
reached at p' = k = \j vn v E c . m . — 0.396 fin . In the insert a discontinuity develops if p' approaches k from 
below or above (see text). 
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FIG. 8: The discontinuity of the real (upper row) and the imaginary (lower row) parts of the off-the- 
energy-shell pure Coulomb t-matrix t c (p,p',x) (given by Eq. (|16p ) near p' = k = \Jm p E c , m , — 0.396 fm -1 
at E l p ab = 13 MeV, p = 0.361 fm -1 , and x = 0.71 as a function of the p' momentum. 
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FIG. 9: (color online) Limiting behaviour of the real part of the off-the-energy-shell screened t-matrix 
t^-(p,p',x) at negative energy E = —13 MeV, p = 0.36 fm -1 , and x — 0.71, as a function of the p 1 
momentum for n = 4 and different values of the screening radius R: R = 20 fm (dotted line) , R = 60 fm 
(dashed-dotted line), R — 120 fm (dashed line), and R — 180 fm (solid line). The pure Coulomb off-shell 
result of Eq. (ri~5]) is given by thick dots. 
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FIG. 10: (color online) The real part of the half-the-energy-shell screened t-matrix t^(p,p, x) before (upper 
row) and after (lower row) the renormalization. The proton lab. energy is E = 13 MeV, and x — 0.71. The 
screening potential was taken with n = 4 and different values of the screening radius R: R = 20 fm (dotted 
line), R = 60 fm (dashed-dotted line), R = 120 fm (dashed line), R = 180 fm (solid line) and R = 500 fm 
(thick solid line) The pure Coulomb half-shell result of Eq. (fT7|) is given by thick dots. 
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FIG. 12: (color online) The real part of the on-the-energy-shell screened t-matrix t^(p,p,x) (upper row) 
and the corresponding renormalized t-matrices (lower row) at = 13 MeV as a function of x for n = 3 
and different values of the screening radius R: R — 20 fm (dotted line), R = 60 fm (dashed-dotted line), 
R = 120 fm (dashed line), R = 180 fm (thin solid line). The Coulomb on-shell amplitude of Eq. (fT8"]) is given 
by thick solid line. The left and the right column differs only in the scale of the x-axis. 
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FIG. 13: (color online) The imaginary part of the on-the-energy-shell screened t-matrix t^}(p,p,x) (upper 
row) and the corresponding renormalized t-matrices (lower row) at E l p ab = 13 MeV as a function of x for 
n = 3 and different values of the screening radius R. For the description of lines see FigfTSl 
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